Equilibrium Shape and Size of Supported Heteroepitaxial Nanoislands 
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We study the equilibrium shape, shape transitions and optimal size of strained heteroepitaxial nanoislands 
with a two-dimensional atomistic model using simply adjustable interatomic pair potentials. We map out the 
global phase diagram as a function of substrate-adsorbate misfit and interaction. This phase diagram reveals 
all the phases corresponding to different well-known growth modes. In particular, for large enough misfits and 
attractive substrate there is a Stranski-Krastanow regime, where nano-sized islands grow on top of wetting films. 
We analyze the various terms contributing to the total island energy in detail, and show how the competition be- 
tween them leads to the optimal shape and size of the islands. Finally, we also develop an analytic interpolation 
formula for the various contributions to the total energy of strained nanoislands. 

PACS numbers: 81.10.Aj, 68.55.Ac, 68.35.Gy 



INTRODUCTION 

The shape and size of adatom islands resulting from growth 
processes has been a subject of numerous recent experimental 
and theoretical studies |EL S, S, 0, S E3, [HI E3, [lH Q [l5, 
I22L 12911 . In the case of Stranski-Krastanow (SK) growth 
mode with islands growing on top of a wetting film, the rel- 
ative abundances of different crystal shapes in a distribution 
show clear peaks at certain volumes, known as optimal sizes 
or magic numbers. In particular, when the growing islands 
are of nanoscopic size the central issue is the possible spon- 
taneous self-organization of islands into arrays with a narrow 
size distribution. Such cases offer immediate technological 
applications in modern nanotechnology. 

In spite of numerous publications in this field published in 
recent years some aspects of such self-assembly process are 
still not quite clear. In particular, there are still uncertainties as 
to whether the observed shapes and sizes of growing islands in 
heteroepitaxy correspond to thermodynamic equilibrium state 
of minimum free energy or limited by kinetic effects. 

In this paper, we examine the minimum energy configura- 
tion of an array of islands with the given constraint of a fixed 
number of adsorbate atoms and with a fixed island density de- 
termined by initial growth conditions J3]. There are various 
physical mechanisms leading to the optimal size and shape of 
the islands. Among them are the relevant surface tensions, 
adsorbate substrate bonding, elastic relaxation in the island, 
wetting film thinning and nonlinear elastic contributions, such 
as bending, buckling and dislocations etc. Ill, 33, 3(1 The 



wetting film and substrate can also mediate elastic dipole in- 
teractions between the islands, which can play an important 
role J3. 

Most of the previous theoretical investigations of this prob- 
lem employ a continuum approach to treat the elastic proper- 
ties of the adsorbate and substrate. In addition, a predefined 
set of facets is usually assumed for the islands together with 



a somewhat arbitrary separation of surface and bulk terms in 
total energy expression. 



The accuracy of these approaches for nanosized islands is 
yet to be determined I37I l39ll . To avoid any such approxima- 
tions, in this work we use a fully atomistic model allowing for 
both elastic and plastic strain relaxation without assumptions 
on predefined shapes. To find the ground state of the atomistic 
systems we apply here the molecular static approach, which 
corresponds to the zero temperature limit. We adopt a two- 
dimensional (2D) model, but extension of many of the results 
presented here to the more realistic 3D systems is also possi- 
ble. The reduced dimension allows to study all the possible 
configurations within feasible computer time. 



This paper is organized as follows. First, we focus on a 
few specific island shapes found for various material parame- 
ters by searching for the minimum energy configuration in the 
small size limit. Then by focusing on these island shapes and 
by varying the misfit and the substrate-adsorbate interactions, 
we find configurations corresponding to the commonly known 
different modes of adsorbate growth, Frank- Van der Merwe 
(FM), Stranski-Krastanow (SK) and Volmer- Weber (VW). For 
the SK mode, we demonstrate the existence of optimal size 
and shape for the islands and study how they vary with the 
coverage of the adsorbate. We discuss the physical mecha- 
nisms that give rise to the narrow distribution of the optimal 
shape and size. Finally, we present an interpolation formula 
for the total energy which yields an accurate description of 
the phase diagram and optimal size and shape of the island in 
comparison with the numerical data. Moreover, this formula 
allows a clear interpretation of the various competing strain 
energy relaxation mechanisms. 



MODEL AND METHODS 



In this work, we adopt a 2D classical atomistic model 
iflii 19 , Eoll . which has been previously used to study mi- 
croscopic mechanisms for strain relaxation in thin films. This 
model also allows for an straightforward extension to the more 
realistic 3D case 12111 . 

Within the model atomistic system is relaxed through stan- 
dard MD cooling procedure and different island shapes are 
systematically compared for finding minimal energy shapes. 
In 2D all relevant island configurations can be studied within 
a reasonable computer time. The model allows for both elastic 
and plastic strain relaxation in heteroepitaxial system without 
any assumptions on predefined shapes of the islands. Some of 
our results here are obtained for relatively small, nanoscopic 
islands up to a few hundred atoms in size, in order to examine 
deviations from the continuum theory of elasticity [Q. 51 0] . 

To allow easy adjustability and anharmonic effects, the in- 
teractions between all atoms in the system are described by 
a modified Lennard- Jones (LJ) pair potential V(r) with 
two parameters, namely the dissociation energy e a b and the 
atomic equilibrium distance r a b, 



Vab(r) = Eab 
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©r ab ,r c (r), (1) 
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r- ab ,r c (O 



1, r< r ab ; 

- 2 (^) 3 ' ^b<r<r c 
0, r c < r. 



(2) 

The modification ensures that the potential and its first deriva- 
tive vanish at the cutoff distance r c . The index (ab) stands for 
substrate-substrate (ss), substrate-film (sf) or film-film (ff) 
respectively. 

The parameter ff s for the adsorbate-substrate interaction is 
simply set as the average of the film and substrate lattice con- 
stants, i.e. rts = (rg + r ss )/2. The lattice misfit / between 
the adsorbate and the substrate can be defined as 



/ = (rff 



OA 



(3) 



A positive mismatch / > corresponds to compressive strain 
and negative / < to tensile strain when the adsorbate island 
is coherent with the substrate. 

As for the interaction parameters, the film-film and 
substrate-substrate interactions are set to be equal, with eg = 
e ss . The film substrate interaction e s f can be parameterized by 
k defined as 



« = (e ss - e sf )/e s 



(4) 




A negative value of k < corresponds to an effectively at- 
tractive and positive k > to a repulsive substrate. 



Figure 1: The generic island shapes A, B, C and D as discovered for 
nanoscopic islands l2(ill . All the shapes have at most a single high 
energy facet. 



Calculations were performed with periodic boundary con- 
ditions for the substrate in the direction parallel to the adsor- 
bate-substrate interface. The size of the unit cell l w (measured 
in units of substrate lattice constant r ss ) determines the den- 
sity of the islands in the SK and VW regime. Two bottom 
layers of the substrate were held fixed to energy minimize a 
semi-infinite substrate while all other layers were free to re- 
lax. In most calculations the thickness of the substrate was 15 
layers, and additional tests were carried out to ensure that the 
thickness did not influence the results. 



RESULTS 

Phase Diagrams in the Submonolayer Regime 

The first issue of interest concerns the equilibrium shapes 
of small islands in the submonolayer regime. In Ref. H20TI 
we presented a complete phase diagram as a function of the 
misfit / and the total number of adsorbate atoms N in this 
regime with k = 0. It was found that there are four charac- 
teristic shapes A, B, C, and D as shown in Fig. Q] Of these 
the shape D, which was not predicted by continuum theory 
calculations 0. was encountered only occasionally for 
larger volumes and misfits. Here, we present the correspond- 
ing results for the cases k = +4% and k = —4%, for N < 85 
and the size of the periodic cell l w = 200. According to our 
tests, the interaction parameter is the most important factor in 
determining the various growth modes in equilibrium. 

Following the procedure explained in Ref. I20I1 . we find the 
ground state shape of an island for a fixed total number Na of 
adsorbate atoms without assuming any predetermined shapes 
by a systematic search approach. Each coherent configuration 
is described by a set of integer numbers, m specifying the 
number of atoms in successive island layers. In terms of these 
numbers, the two types of facets, considered in the previous 
works 10, 13] correspond to rij — nj+i = 1 for steep facets and 
rii — ni + i = 3 for shallow facets. The only physical restric- 
tions we impose are that the island has a reflection symmetry 
about a line through the center and that overhangs are not al- 
lowed. Then, for each initial configuration, molecular dynam- 
ics (MD) cooling is run to allow the system to relax and reach 
a minimum energy configuration. The equilibrium shape for a 
given Na is identified as the relaxed island configuration with 
lowest energy among all the configurations. 
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Figure 2: Phase diagrams for submonolayer regime (Na < 85) with 
l w = 200. Panel (a) is for k — +4% and (b) for the case n — 
—4%. The line in panel (a) is a small domain of configurations that 
belong to the Frank- van der Merve regime (FM). As expected, in the 
attractive substrate case island formation and re-entrant transitions 
are suppressed for fixed / when Na increases. See text for details. 



The phase diagram for the case where the substrate is re- 
pulsive with k = +4% is shown in Fig. |2|a). As expected, 
island formation is enhanced as compared to the case where 
K = 0. Even with zero misfit, island formation starts around 
Na ~ 50 and the wetting film disappears completely beyond 
Na ~ 75. We detect no new island shapes or any occurrences 
of the shape D. 

The results above should be contrasted to the case of an 
attractive substrate with k = —4% in Fig. Hb). Complete 
wetting regime dominates the phase diagram almost up to the 
largest values of mismatch considered here. As for the shape 
of the islands, all but one case belong to the categories A and 
C. The single exception at / = —10% and Na = 80 falls in 
the category D. We also note that as in Ref. 1201. the phase 
diagrams here are asymmetric with respect to tensile and com- 
pressive strain highlighting the importance of the anharmonic- 
ity of the atomic interaction potentials. 



Global Phase Diagram at Higher Coverage 

At higher coverages with N 3> l w , it is possible to have 
in addition to the partial wetting film (FM) and small island 
(VW) phases also the phase where islands grow on top of wet- 
ting films (SK). 

We have calculated the global phase diagram (GPD) as 
shown in Fig. [3]for Na = 820 and l w = 200 (coverage ~ 4 
layers) in the k — f plane, where |/| < 7% and |k| < 6%. 
Numerical data for this diagram have been taken with steps of 
0.5% along both axes. 

Following the approach of Scheffler et al. SHU 3, the 
growth process can be divided into different stages. The early 
nucleation stage mainly determines the island density. In our 
model, this is fixed by the size of the unit cell l w which is 
inversely proportional to the island density. The second stage 
is where the islands can grow at the expense of the wetting 
layers. The main driving force here is a pathway to lower the 
strain energy while not sacrificing too much the adsorbate- 
substrate interface bonding energy. 

The data in Fig. [3] are calculated by the following steps. 
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Figure 3: Global phase diagram of island shapes as a function of the 
misfit parameter / and the interaction parameter k for TV = 820. For 
each point in the diagram the minimization was done with fixed k 
and /. The three different phases corresponding to the well known 
growth modes are indicated in the diagram. See text for details. 



For a fixed Na at 820, we investigated all configurations cor- 
responding to an island with shape corresponding to A, B or 
C as shown in Fig. [TJ together with any number of complete 
wetting film layers and at most one partial wetting film layer 
(layers which are more than seven atoms wider than the is- 
land base are counted as partially filled wetting layers). The 
initial configuration has the equilibrium lattice constant of the 
substrate. When the energy of the system is minimized with 
molecular dynamics cooling, the adsorbate releases strain en- 
ergy by relaxation. This local minimization technique does 
not lead to plastic deformations in this parameter range. Each 
point of Fig. [^corresponds to a relaxed configuration with the 
lowest energy. 

As can be seen in Fig. [3j two well-known phases already 
identified in the submonolayer regime are present in different 
regions of the phase diagram for the higher coverage situation: 
the Frank-van der Merwe (FM) phase dominates for relatively 
small misfits and attractive substrate interactions while the 
Volmer- Weber (VW) phase is favored for larger misfits and 
posive values of n. As the lattice misfit / increases, there is a 
transition from wide, relatively thin islands to sharper forms. 

However, perhaps the most interesting feature of the phase 
diagram for coverages beyond two layers is the appearance 
of the Stranski-Krastanov (SK) regime for large misfits and 
negative values of k. In the SK regime, islands exist on top 
of wetting films. With the present interaction parameters, 
the SK phase lies between |«| > 1% and / > 4%. The 
triple points separating the phases are located approximately 
at (kJ) = (0.5,4)% and («,/) = (-0.5,-3.5)%. Again, 
the phase diagram is asymmetric between tensile and com- 
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Figure 4: A schematic figure of a typical island geometry considered 
here. The height of the island and the wetting film are denoted by hi 
and h w , respectively. The island-island separation, which equals the 
width of the periodic simulation cell is denoted by l w , and the island 
base width is U. 



pressive strain values, even for relatively small values of /. 

We have numerically calculated the global phase diagrams 
for other values of Na — 300 and 1000 and the results are 
very similar. For coverage Na = 300, the phase boundary 
between the FM and SK regimes becomes less well defined. 
However, for Na = 1000, the changes to Fig. [3]are minor; the 
positions of the triple points seem to move to smaller values 
of l/I- 



Optimal Shape and Size for Islands in the SK Mode 

For technological applications, it would be desirable to 
have a self-organized array of islands which has a very narrow 
size and shape distribution function. Experimentally, such 
"optimal" island sizes have been observed in some heteroepi- 
taxial systems 113 ll 13211341 13511 . The basic question is whether 
the optimal size islands correspond to equilibrium minimum 
energy configurations, or are they just a consequence of ki- 
netic limitations that prevent further growth on accessible time 
scales. The existence of optimal shapes from purely energetic 
considerations has been prov en for some specific systems such 
as those in Refs. UmlHIIil. In Fig. |3 we define a num- 
ber of parameters which characterize the shape and size of 
the island and the wetting film. The height of the island and 
the wetting film are denoted by hi and h w , respectively. The 
island-island separation, which equals the width of the pe- 
riodic simulation cell is denoted by l Wi and the island base 
width is ij. A schematic diagram is shown in Fig. [4] 

We have investigated the optimal shapes and sizes of the is- 
lands in the SK mode for the misfit parameter / = 5% and 
interaction parameter k = —4% as a function of different 
amounts of adsorbate JV^, for an island density correspond- 
ing to l w = 200. The results are shown in Fig. [5] At low 
coverages 8 = Na/I w the atoms completely wet the substrate 
and form wide, non-wetting islands after which the actual SK 
regime emerges. For coverage 6 > 4, both the base width 
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Figure 5: Figure demonstrating the appearance of optimal island 
shapes and sizes in the SK regime. The direct energy minimization 
data points, which are marked with dots, are obtained by minimizing 
the energies of all acceptable shapes of a given coverage and select- 
ing the lowest energy configuration among them. The solid line is an 
estimate given by numerical minimization of the interpolation for- 
mula. 



li and the aspect ratio hi/U of the island increase slowly with 
coverage. Typically, the adsorbate forms a combination of one 
full and one incomplete wetting layer together with an island 
of shape A. 

To understand the strain relaxation in the SK phase, we se- 
lect the 9 = 2.5 point from Fig. [5] corresponding to total 
number of adatoms Na = 500. The optimal shape and size 
of the island in this case has li = 41 and hi = 8 on top of 
one complete layer of wetting film, see panel (b) of Fig. [6] 
In addition to this configuration, we show the spatial energy 
distribution of the totally relaxed configurations for two other 
geometries selected from the high and low aspect ratio ends 
of the shape spectrum, namely configurations with Z, = 24, 62 
and hi = 24, 5, respectively. The strain patterns of these dif- 
ferent configurations are shown in panels (a) and (c) of Fig. 
[6] We note from these three geometries that the strain pattern 
falls into two different categories. For the large aspect ratio 
(li = 24, hi = 24), except for the immediate adsorbate sub- 
strate interface region, the bulk of the island is largely relaxed. 
In the other limit of a small aspect ratio (Z, = 62, hi = 5), ex- 
cept for the surface and corner region, most of the atoms inside 
the island are strained. For small islands, the relaxation fa- 
vors sharp islands but the adsorbate-substrate interface bond- 
ing energy and the surface energy favor wide islands. It turns 
out that the most efficient way of lowering the strain energy 
while still maintaining a strong interface bonding energy is to 
pick a geometry in between these two extremes. This balance 
determines the aspect ratio and the optimal shape of the is- 
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Figure 6: The qualitative change in the relaxation pattern in wide and 
narrow islands. These shapes are taken from a larger configuration 
set which was used to fit Eq. l[5}. The darker colors indicate in- 
creasing levels of strain. To enhance contrast, the outmost adsorbate- 
vacuum surface layer is not shown. The configurations shown here 
correspond to Na — 500. The islands have shape A on top of one 
complete layer of wetting film. The widths are 24, 41 and 62 from 
top to bottom and the heights hi are 24, 8 and 5, respectively. The 
configuration in the middle with aspect ratio hi/U = 0.2 and base 
width U = 41 has the lowest energy. 



land. When the coverage increases, the total strain energy can 
be lowered by allowing the size of the island to increase at the 
expense of reducing the thickness of the wetting film. The op- 
timal size of the island results from the balance between the 
strain relaxation of larger island and the "thinning" energies 
of the wetting film. Another factor that contributes to deter- 
mining the optimal size of the island is the indirect interac- 
tion between the islands mediated by the substrate relaxation. 
These points are discussed in more detail in next section. 



ANALYTIC INTERPOLATION FORMULA FOR TOTAL 
ENERGY 

In this section, we will look in more detail into the various 
competing mechanisms leading to the existence of an opti- 
mal size and shape for the islands in the SK regime. We also 
develop an approximate analytic interpolation formula that in- 
corporates all these mechanisms and allows a simple interpre- 
tation of our numerical data. Ou r app roach is in many respects 
similar to those in Refs. ifHil^lzl |27ll . 



To develop an analytic interpolation formula, we first 
choose a reference system corresponding to a completely un- 
strained substrate and completely strained adsorbate such that 
the in-plane spacing between atoms in each layer of the island 
is the same as that in the substrate. This reference system has 
three kinds of bonds: unstrained, strained and substrate-ad- 
sorbate interface bonds. The corresponding energies are e u b 
for unstrained bonds and e xa and e ya for horizontal and ver- 
tical adsorbate bonds, respectively. The substrate-adsorbate 
interface bonds, which are vertical, are denoted by e y i. These 
are obtained from Eq. ((TJ, and the details are given in the 
Appendix. 

The bulk contributions Eb and Ew for the substrate and 
adsorbate are proportional to the number of atoms in the ad- 
sorbate and substrate, respectively. In addition, each surface 
and interface leads to a correction term due to missing or 
changed bond energies, resulting in a total surface energy Eg. 
For the wetting film, partial layers are allowed to account for 
the situation where the number of atoms in the wetting film de- 
viates from multiples of l w . All corner energies are neglected. 
The different contributions to Eb , Ew , and E$ in terms of the 
bond energies e u b, e xa , e ya , and e yi are listed in Table Q] above. 



Relaxation Energy of Island 

The total energy of a fully strain relaxed island can now 
be written as E'totai = E xe { + f2, where fl is the relaxation 
energy. As the panel (c) of Fig. (O shows, the relaxation is 
seen only in the vicinity of the non-horizontal facets when the 
island aspect ratio is low. In this situation the top facet and the 
interior are both completely stressed. When the qualitative 
pattern of relaxation is of this kind, the energy difference £1 
has to become practically independent of the island width for a 
set of islands with the same height. If we consider the islands 
with constant height in the high aspect ratio case, the panel 
(a) of Fig. (|6]l shows that the relaxation energy should depend 



Table I: Reference energy expressions. The quantity lo — l w — h 
when h w > 1 and lo — If when h w < 1 and If is the width of 
the partial wetting fraction, which is not under the island. The sym- 
bols e u b , e xa and e ya refer to unstrained substrate bond and strained 
horizontal and vertical adsorbate bond energies, respectively. The 
symbol e y i is the energy of a strained substrate-adsorbate bond. 



Contribution 
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Substrate bulk Eb 


— 6e u biVs 


Adsorbate bulk Ew 


-(2e xa + 4e ya )iVA 


Surfaces Es : 




Substrate- vacuum 


2e u b(lw — li — lo) 


Wetting film- vacuum 


2e ya 'o 


Oblique facet- vacuum 


2(^ y a ~t~ ^xa)^i 


Island top-vacuum 


2c y a (li h% ) 


Substrate-adsorbate 


(2e ub + 2e ya - 4e yi )(fo + li) 
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almost linearly on the width. 

To quantitatively verify this, we study a set of islands with 
different values of island base widths ^ and island heights hi. 
The wetting layer thickness h w and the ratio l w /U are kept 
constant with h w = 2 and l w /U = 10. This value of l w /k 
is large enough so that the interaction energy between islands 
is negligible at this separation. The chosen island heights hi 
are 5, 15, 31, 45 and 59. The base widths U for each height hi 
are chosen with a denser mesh near the large aspect ratio limit 
and with a sparser mesh in the small aspect ratio, shallow is- 
land limit. The energy of the island in each case is minimized 
and then the additional relaxation energy is obtained after the 
subtraction of the reference energy. The numerical data of the 
extra relaxation energy for the entire set of islands with differ- 
ent aspect ratios can be fitted to a simple analytic form which 
has the properties described above, namely 

« = ^shallow -Kk-lc- V(ih-lc) 2 +d 2 ) . (5) 

The first term Shallow in Eq.(5) represents the relaxation en- 
ergy in the shallow island (small aspect ratio) limit, which can 
be fitted to a form independent of lj given by 

-^shallow -0.194(7)^(^ + 8.16(7))^. (6) 

Here the term in hf is related to the volume relaxation and 
the term linear in hi comes from the surface relaxation of the 
shallow island. It should be noted that when the island base U 
is in the regime where depends only weakly on li and the 
height hi becomes large, Eq. © breaks down. In this large 
volume limit the term is linearly dependent on hi, i.e. it has 
turned into a correction to the surface energy of the vertical 
island facets. 

The other variables l c , b, and d of Eq. (f5]) are functions 
of the height of the island and the numerical fit of these vari- 
ables are given in Table [TT1 The value of l c defines a crossover 
between the shallow island limit and the sharp island limit 
where the qualitative nature of the relaxation pattern in the is- 
land changes. For a fixed hi, the energy of the island starts to 
increase sharply when li is reduced below the threshold value 
l c which is a function of hi. 

To keep the analysis simple, we have not included small 
contributions to the total energy, such as corner energies and 
stretching energy of the adsorbate-substrate interface in the 
interpolation formula. There exist rare circumstances when 
these effects can become important and in these cases the for- 
mula does not hold. 

It should be noted that at some coverage the island volume 
becomes necessarily so large that the relaxed neighborhoods 
of the left and right facets do not overlap even in the equilat- 
eral case. In this limit the relaxation energy is proportional 
to the length of the non-horizontal facets and can be adsorbed 
to the corresponding surface tensions. In this case the above 
formula is not expected to hold. 

A similar study concerning the stability of strained het- 
eroepitaxial systems has been conducted in Ref. 12711 . They 
find that an island array is unstable in all the cases considered. 



Table II: Parameters of Eq. (f5]( as function of height. 



Parameter 


Expression 
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0.0134(2) (/m + 1.602(2)) 


lc 


2.45(7) A< + 8.47(3) 


d 


1.21(0)fti +8.82(0) 



The main difference to our study is that we do not assume the 
island-island distance to be fixed by the interaction between 
the islands but by the nucleation stage, as in Refs. 0,0,0]. 

Comparison of Numerical Data with Interpolation Formula 

In this section, we will use the interpolation formula ob- 
tained in the previous section to examine the global phase dia- 
gram and determine the optimal sizes and shapes of the islands 
in the SK regime. First, we examine the global phase dia- 
gram. For this purpose, we fix the total number of adsorbate 
atoms Na = 820 and determine the minimum energy config- 
uration as a function of the misfit parameter / and interaction 
parameter k by comparing the total energies of different con- 
figurations. Note that the reference energy E ie { already has 
a dependence on k and /. For the additional relaxation en- 
ergy O we take the formula, which was fitted for / = 5%, 
K = —4%. We verified numerically that in the first approx- 
imation it scales with the misfit as / 2 and neglect its depen- 
dence on k Roll . The resulting GPD as shown in Fig. |7]is very 
close to the one obtained from the direct energy minimization 
as shown in Fig [3] The disagreements have generally three 
sources. 

First, the interpolation formula is constrained to islands of 
shape A in some points of the phase diagram. This can lead 
to two kinds of differences. As we showed earlier, when 
all configurations are compared with the direct minimization, 
the lowest energy shape is almost always among A, B, or C 
shapes so that the extra facets in the B and C shapes are of 
the shortest possible length. However, at some coverages the 
combinatorics limits the number of these three shapes to be 
fairly small. It can happen that these shapes are also of ex- 
tremal widths or thicknesses whose leading energy contribu- 
tions are out of the average range. Under these circumstances 
the configuration space would be sampled better by allow- 
ing the shapes to have some asymmetry or additional facets 
at the corners. This finite size effect is relevant for small is- 
land volumes only, and as we showed it is not very signifi- 
cant even in this regime. Anyhow, if the interpolation formula 
is minimized analytically, the configuration space is continu- 
ous and does not suffer from this finite size effect. Another 
kind of difference originating from the same source is the fact 
that sometimes the small terms (such as the corner energy) 
not present in the formula can force the direct energy mini- 
mization to switch to a wider island width at a lower coverage 
than what the interpolation formula gives. This can reduce 
the island height by one and increase the width by two layers 
but the effect is not cumulative and cannot separate the direct 
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minimization and interpolation results far apart. 

Second, close to the critical coverage separating two growth 
modes the energies of very different configurations, such as a 
flat layer and an island, are comparable. In this kind of situ- 
ation the small energy terms not present in the interpolation 
formula can actually determine the energy balance and the 
real lowest energy shape is therefore not detected by the for- 
mula. However, when the coverage is outside the transition 
lines these cases become irrelevant. 

Third, the analytic formula is less reliable in the large aspect 
ratio limit (sharp islands) with high coverages. The height de- 
pendence of the parameters is fitted for the small island limit 
where the surface and bulk energetics are not clearly sepa- 
rable. With sufficiently large island volume the relaxation 
becomes clearly associated with the island facets and in this 
case the relaxation energy must be directly proportional to the 
length of the associated facets. Since in this case the stress 
energy cost increases proportionally to the island volume, dis- 
locations nucleation becomes relevant; an effect which is not 
included in the present model. In this limit also the validity of 
the f 2 scaling of the relaxation energy breaks down. 

Finally, we should also mention that the island-island in- 
teraction energy is neglected in the present interpolation for- 
mula. However, we have checked that this effect is negligible 
at a unit cell length of l w = 200 and does not affect the results. 

Next, we investigate the optimal size and shape of the is- 
lands in the SK regime for the parameters / = 4% and 
k = —3%. For this purpose, we take the interpolation for- 
mula for O derived for a fixed thickness for the wetting film 
(h w = 2) and apply it to configurations for other values of 
wetting film thickness. This amounts to neglecting the wet- 
ting film thinning energy and is only expected to be valid in 
the low coverage limit. We vary the size and aspect ratio of 
the island for a given coverage through choosing different val- 
ues for li and hi- The corresponding value for the thickness of 
the wetting film h w is then determined by the constraint of the 
total amount of the adsorbate material being constant. The to- 
tal energy expression is minimized to yield the optimal value 
of k and hi, which determine the optimal size and shape for a 
given coverage. 

In Fig. [5] we show the island base width li and aspect ra- 
tio hi/li as a function of the coverage. The predicted optimal 
island base widths from the direct numerical calculation and 
those obtained from the total energy minimization with the in- 
terpolation formula are close. The optimal aspect ratios deter- 
mined from the direct minimization and from the interpolation 
formula have the same qualitative coverage dependence. 



SUMMARY AND CONCLUSIONS 

In this work, we have investigated the various phases re- 
sulting from the relaxation of strain energy starting from ini- 
tial epitaxial strained adsorbate wetting films. We employ the 
molecular static method which is a powerful tool to overcome 
the restrictions related to conventional minimization methods 
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Figure 7: The analytic interpolation formula for the total energy 
yields a global phase diagram in good agreement with one obtained 
from direct numerical minimization. The fitting of the formula was 
done in the point where / = 5% and k = —4%. The small asymme- 
try of the direct minimization is related to the asymmetry of the inter- 
action potential which is especially pronounced near the locations of 
the largest displacements such as the corners, as Fig.[6]demonstrates. 
In the large scale the direct minimization result is fairly symmetric, 
which suggests that the rest of the relaxation energy scales roughly as 
f 2 as the continuum theory predicts. Since the interpolation formula 
already neglects the corner energies, we also take the elastic energy 
scaling to be of the f 2 form. This explains the mirror symmetry of 
the positive and negative misfit sides of the interpolated phase dia- 
gram in the lower panel. 



that rely on the continuity and differentiability of the underly- 
ing manifold. In this work the phase diagram of 2D heteroepi- 
taxial thin film systems was obtained by direct minimization 
for various combinations of the lattice misfit / and parame- 
ter k describing the interface tension. In the submonolayer 
regime, the Stranski-Krastanow (SK) phase is not present, but 
a transition between layer-by-layer (FM) equilibrium shapes 
and Volmer- Weber (VW) shapes is sensitive to the value of 
k. As the adsorbate coverage increases, the SK phase appears 
between the FM and VW regions when k had a value that 
was favorable for the contact between substrate and adsorbate. 
The islands always have truncated pyramid shapes with minor 
edge defects. 

In the SK regime, we have investigated the important issue 
concerning the optimal shape and size of islands by searching 
through the minimum energy configurations of systems for 
fixed material parameters / and k. The shapes were restricted 
to truncated pyramids. The actual optimal size and aspect ra- 
tio depend on the coverage of the adsorbate atoms. Both the 
aspect ratio and the base width of the island increases slowly 
with coverage. To understand better the numerical results, an 
analytic interpolation formula was developed. This approxi- 
mate formula takes into account the strain relaxation from a 
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reference state of totally strained wetting films, adsorbate is- 
land and an unstrained substrate. It produces results that agree 
very well with the numerical energy minimization data. 

The optimal shape and the aspect ratio result from a com- 
promise of increasing the relaxation energy for a sharp (higher 
aspect ratio) island vs. increasing the adsorbate-substrate in- 
teraction for a shallow (low aspect ratio) island. The optimal 
size of the island grows as the coverage increases because a 
larger island is more effective in strain relaxation and lower- 
ing of the elastic energy. Eventually, when the base width of 
the island ij grows to be comparable with the unit cell size l w , 
the optimal island size will be determined by the substrate me- 
diated indirect interaction between the islands, as suggested 
in Ref. Il24ll . It is important to note that the details of the 
coverage dependence of the optimal size in our model differ 
somewhat from the results of Ref. 17[], where the wetting film 
thinning energy was evaluated using a continuum approach 
and is finite at thickness even above two layers. In our model, 
the wetting film thinning energy implicitly present in both the 
interpolation formula and the numerical energy minimization 
is negligible until the last layer of wetting film starts to get de- 
pleted. Thus the optimal size at low coverages shown in Fig. [5] 
corresponds to the total adsorbate atoms minus one complete 
plus possible one partial wetting film layer. 

Finally, we would like to discuss the applicability of 
our results for realistic 3D cases. While the fact that our 
global phase diagram correctly reproduces all the well-known 
growth modes in the relevant regimes indicates the overall va- 
lidity of the model, there are several additional features in 3D, 
which need to be considered. Assuming no intermixing of the 
substrate and film atoms, the most important ingredient miss- 
ing in 2D is the spatial anisotropy of various crystal surfaces. 
Recent experiments in 3D indicate II34I1 that changing the mis- 
match of the substrate can induce significant morphological 
changes in 3D islands. For the close-packed (111) surface ge- 
ometry, which most closely corresponds to the 2D case stud- 
ied here, we have done minimization of selected island shapes 
and find that there are stable island configurations, which are 
straightforward 3D generalizations of the 2D shapes found 
here. However, a systematic search for optimal 3D nanois- 
land shapes is computationally very expensive and beyond the 
scope of the present work. 
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APPENDIX: BOND ENERGIES IN EPITAXIAL STRAINED 
FILM 

An adsorbate bulk atom has two horizontal and four diago- 
nal bonds, the energies of which are denoted by e xa and e ya , 
respectively. The horizontal distance of the atoms in a co- 
herent adsorbate layer is r ss and the equilibrium distance is 
x' = r ss (l + /)■ The energy of a horizontal bond — e xa in the 
reference state is thus Vff(r ss )/2. In a hexagonal lattice with 
pair potential the Poisson ratio is always v = 1/3. Thus, the 
ratio of vertical and horizontal strains u ml and u TT is 



^yy 



or 



V_ 

v' 



3x' n 



(7) 



(8) 



where y' is the actual distance between adjacent horizontal 
layers in a strained adsorbate and y' is the corresponding equi- 
librium distance, y' = \/3x' /2. This gives 



, _ , / 4x' r s 
V ~ Vo \ 3x' 



(9) 



The length of the diagonal bond in the initial strained film 
is given by the relation r 2 d = r"l s /A + 3(y') 2 /4 and energy of 
the bond — e ya is given by Vg/2 of Eq. (Q]) with the argument 
rd- 

The bond strength between the substrate and adsorbate e y i 
can be easily obtained from the adsorbate bond e ya with two 
modifications. First, the misfit parameter / is replaced by f/2 
and then the overall strength is multiplied by a factor 1 + k. 
The first modification arises from the fact that the equilibrium 
distance between the substrate and adsorbate atom is the arith- 
metic mean of the equilibrium distances between pure mate- 
rials and the second follows from the definition of k relating 
the strength of the film-substrate interaction to the film-film 
interaction. 
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